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CONVERGENCE OF A FULLY DISCRETE VARIATIONAL SCHEME FOR 

A THIN-FILM EQUATION 


HORST OSBERGER AND DANIEL MATTHES 


Abstract. This paper is concerned with a rigorous convergence analysis of a fully discrete 
Lagrangian scheme for the Hele-Shaw flow, which is the fourth order thin-fllm equation with 
linear mobility in one space dimension. The discretization is based on the equation’s gradient 
flow structure in the L^-Wasserstein metric. Apart from its Lagrangian character — which 
guarantees positivity and mass conservation — the main feature of our discretization is that it 
dissipates both the Dirichlet energy and the logarithmic entropy. The interplay between these 
two dissipations paves the way to proving convergence of the discrete approximations to a 
weak solution in the discrete-to-continuous limit. Thanks to the time-implicit character of the 
scheme, no CFL-type condition is needed. Numerical experiments illustrate the practicability 
of the scheme. 


1. Introduction 

1.1. The equation and its properties. In this paper, we propose and study a fully discrete 
numerical scheme for the following nonlinear fourth order equation on the bounded domain 
n = (a, 5), a < 6 with no-flux boundary conditions: 

dtu = — (uuxxx)^ + (Vxu)^ for t > 0 and x e O,. (1) 

uux = 0, uuxxx = 0 for t > 0 and x e dO. (2) 

We assume that the potential V e C‘^{Q) is non-negative with bounded second derivative, 

V ^ 0, A := sup \Vxx\ < 00 , (3) 

a typical choice being V{x) = ^x'^. Equation 0 belongs to the family of thin film equations, 

dtU = — div(m(u) D Au). (4) 

Equations of this form are used as reduced models for laminar flow with a free liquid-air interface 
m- The case of linear mobility m{u) = u studied here is further used to describe the pinching 
of thin necks in a Hele-Shaw cell, hence Q is often referred to as the Hele-Shaw flow. 

The analysis of the fourth order degenerate problems in Q is by no means trivial. There 
exists a rich literature on the existence and long-time behavior of solutions, see e.g. [iimiiiiii 
[3 uni HD- These results could not be obtained by straight-forward extension of the techniques 
for second order parabolic equations. The most decisive difference between second and fourth 
order is the lack of comparison principles for the latter. Instead, energy and entropy methods 
play a key role in the analysis. Naturally, classical parabolic theory applies in zones on which the 
solution u is strictly positive, yielding C”-regularity there. However, one is typically interested 
in solutions that are not strictly positive but have a compact, time-dependent support. For such 
solutions, one only has the regularity that is induced by the energy/entropy estimates, which is 
usually something of the type L®([0, T]; n L^([0, T]; R^(r2)), but no better. 

It is known that Q carries a variety of structural properties: solutions conserve mass and 
non-negativity, there exists a large class of Lyapunov functionals |25j , and it can be written as a 
gradient flow in the energy landscape of the following (modified) Dirichlet functional, 

8^(u) = £(u)V(u) , with £{u) = - f (^dxu)^ dx, V{u) = ( V{x)u{x)dx, (5) 

2 Jq Jn 

with respect to the L^-Wasserstein metric [15] • The two main a priori estimates for the well- 
posedness theory of Q are obtained from the dissipation of £^, and from the dissipation of an 
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auxiliary Lyapunov functional, namely the entropy, 


T-l}^^{u) = / uloguda;. 

Ja 


Formally, the respective estimates are immeditaly obtained from an integration by parts; for 
V = 0, one obtains 


= [ UxiuUooxx)xxdx 

Jn 

d f 

= / ^Ogu{uUxxx)xdx 

dt Jn 


= / uulxxdx, ( 6 ) 

Jn 

= I ulxdx. ( 7 ) 

In Jn 

Notice that energy dissipation does not provide L^([0, T]; i7^(n))-regularity, due to the degen¬ 
eracy of the integrand in regions where u vanishes. In principle, the famous Bernis estimates [3] 
could be used to extract an L^([0, T]; i7^(n))-bound on but we shall not discuss this ansatz 
here since these — quite delicate — estimates seem impossible to discretize. 

In the numerical approximation of solutions to 0 , it is natural to ask for a structure¬ 
preserving discretization that inherits at least some of the properties listed above. A minimal 
criterion is to guarantee nonnegativity of solutions — which turns out to be already a challenge. 
Here, we try to push the structure preservation as far as possible, with a scheme that translates 
both the Lagrangian and the gradient flow nature of Q from continuous to discrete, and which 
inherits simultaneously the Dirichlet energy and the entropy as Lyapunov functionals. We even 
obtain a discretized version of the estimate 0. and that is the key element for the convergence 
analysis. 

Our discretization ansatz is closely related to the one that has been introduced and analyzed 
recently by the authors [29j in the context of the fourth order quantum drift diffusion (QDD) 
equation, also known as Derrida-Lebowitz-Speer-Spohn (DLSS) equation: 

dtu = -{u{\ogu)xx) + {yxu)^. ( 8 ) 

This equation is a gradient flow in the L^-Wasserstein metric as well (TH]. In place of , the 
flow potential is given by the perturbed Fisher information 

= JF{u) + V{u), with JF{u) = 2 / [dxVuY dx, 

Jn 

and V as above. There exists a non-obvious connection between ([^ and the linear heat equation 
m, that is best understood as a relation between iF, and the L^-Wasserstein metric p7] . 
The main feature of the particular discretization of (|^ used in [55] is the preservation of that 
relation, and that paves the way to a relatively complete analytical treatment of the resulting 
numerical scheme. In the context at hand, the connection between £ and — that is, the 
origin of the relation Q — is less obvious, but on grounds of the ideas developed in [29] , we are 
able to define a special discretization that admits a discretized version of Q. 

1.2. Definition of the discretization. We are now going to present a discretization for Q 
that approximates weak solutions to Q of finite positive mass M e K>o. Basically, we follow 
the ansatz from [2^, but we shall deviate in the discretization of the potential of the flow. First, 
the equation Q is re-written in terms of Lagrangian coordinates: since each u{t,-) is of fixed 
mass M, we can introduce time-dependent Lagrangian maps X(t, •) : [0, M] —> H implicitly by 

rnt,0 

^ = / u{t,x)dx, for each ^ e [0, M]. (9) 

J a 

For the moment, we ignore the ambiguity in the definition of X(t,^) outside of the support of 
u{t). Expressed in terms of X, and after elementary manipulations, the Hele-Shaw equation Q 
becomes: 


d,X = d,[-Z^ 


+ 4^' 


-t- F(X), where Z(t,C) := = w(t,X(t,0). 

^ ’ (10) 
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It is easily seen that equation (10) is the L^-gradient flow for 

2 


E 


1 

(X) ■.= £^{uoX)=- I 


1 


X, 




nM 


V{X)d^, 


with respect to the usual L^-norm on L^([0, M]; R). This directly reflects the gradient flow 
structure of ( [I|) w ith respect to the L^-Wasserstein metric. 

Equation (10) is now discretized as follows. First, fix a spatio-temporal discretization pa¬ 
rameter A = where r > 0 is a time step size, and S = M/K for some K e N defines 

an equidistant partition of [0,M] into K intervals of length 6 each, i.e., = kS for 

k = 0,1,..., K. Accordingly, introduce the central first and second order finite difference op¬ 
erators DJ and for discrete functions defined either on the ^fc’s or on the interval midpoints 
^fc+i /2 = {k + 1/2)5 in the canonical way; see Section 


2.1 


for details. 


At each time t = nr, the Lagrangian map X{t, •) is approximated by a monotone vector 
x2^ = {xi, ..., x'^_i) e with o < x” < • • • < x^_i < h 

in the sense that X{nT, kS) x^. We will further use the convention that xq = a and xk = b. 
For brevity, introduce the vectors = z[x2^] with entries 

<5 1 


- 


[Dji], 


for K = 


1 3 


2 ’ 2 ’ 




( 11 ) 


and z_i = zi and = Zj(_i by convention. These vectors approximate the function Z in 

(101 such that Z{nT,KS) ^ z”. 

The fully discrete evolution for the x^ is now obtained from the following standard discretiza¬ 
tion of (10) with central finite differences: 

„n—1 

'(z-)2D2[(g-)2] +VM- (12) 




= DJ 


-(z")3 
2^ ’ 


D?rz" 


1 


Note that there are infinitely many equivalent ways to re-write the right-hand side of equation 
(101, and accordingly infinitely many (non-equivalent!) central finite-difference discretizations. 
Another one, having different properties, is studied in [31] . Our convergence result only applies 
to the particular form (12), since only for that one, we obtain “the right” Lyapunov functionals 
that provide the a priori estimates for the discrete-to-continuous limit. 

Finally, we define a time-dependent, spatially piecewise constant density function : O —» 
RjtO from the sequence xa := (^)”=o via 

ul = U5[x^] := ^ (13) 

By definition, the densities are non-negative and of time-independent mass, 

p nXp. c 

/ uldx=Y, / = K5 = M. 


In 


k = l • 


Finally, we introduce the piecewise constant interpolation {ma},- ^ KlgtO ^ H 

{uaIt. (t) = u’a (^ ~ 1)'^ < t ^ nr, 

and {maIt. (0) = u%. 


in time by 
(14) 


1.3. Main results. For the statement of our first result, fix a discretization parameter A = 
(r; 5). On monotone vectors x e R^“^ with densities z = z[x], introduce the functionals 

K r .ff \ 2 K 


H5^^(x) log(^fc-i). 

k=l 


Er 


(x) := „ 2 


Zuii H” Zu_i / Zh.\i — Zh_i 

rv Q M. Q I Ay Q Ay Q 


fc = l 


+ <5 2 K(Xfe), 


fc =0 


which are discrete replacements for the entropy and the modified Dirichlet energy functionals, 
respectively. 
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Theorem 1. From any monotone discrete initial datum a sequence of monotone satis 
fying (12) can be constructed by inductively defining ^ as a global minimizer of 


^)2 + E ^( x ). 


(15) 


This sequence of vectors dissipates both the Boltzmann entropy and the discrete Dirichlet 
energy, 

and £^(5^) E,^(x^-i). 


To state our main result about convergence, recall the definition (14) of the time interpolant. 
Further, A symbolizes a whole sequence of mesh parameters from now on, and we write A ^ 0 
to indicate that r —» 0 and (5^0 simultaneously. 

Theorem 2. Let a non-negative initial condition e H^(n) of finite second moment be given 
and fix a time horizont T > 0. Choose initial approximations x^ such that = u[x^] ^ 
weakly in as A ^ 0, and 


W := supE^(j4) < 00 , Htd := supHl5^^(34) < oo. 

A A 


(16) 


For each A, construct a discrete approximation xa according to the procedure described in Theo- 
rem^above. Then, there are a subsequence with A —> 0 and a limit function e C'([0,r] x S2) 
such that: 

• {'ua}t converges to locally uniformly on [0,r] x fl, 

• w*(0) = ■a°, 

• M* satisfies the following weak formulation of ([^ with no-flux boundary conditions ([^.' 


/ / dtppu,^ dtdx + / riN{u,^,p)dt = 0, 

Jo Jr Jo 


(17) 


for any test functions p e C'®(n) with p'{a) = p'{0) = 0, and p e C'®((0,T)), where the 
operator N is given by 

N{u,p) ■=7;[ {{ir)xPxxx + SuIpxx) dx+ f V:rUPxdx. (18) 

^ Jn Jn 

Remark 3. (1) Quality of convergence: Since {ma},- is piecewise constant in space and 

time, uniform convergence is obviously the best kind of convergence that can be achieved. 

(2) Rate of convergence: Numerical experiments with smooth initial data show that the 
rate of convergence if of order t + 5^, see Section 

(3) No uniqueness: Since our notion of solution is very weak, we cannot exclude that different 
subsequences o/Ima},- converge to different limits. 

(4) Initial approximation: The assumptions in ( 161 are not independent: boundedness of 
E^(x^) implies boundedness o/H^^^(x^) from above. 


1.4. Relation to the literature. The idea to derive numerical discretizations for solution of 
Wasserstein gradient flows from the Lagrangian representation is not new in the literature. A 
very general (but rather theoretical) treatise was given in [53]. Several practical schemes have 
been developed on grounds of the Lagrangian representation for this class of evolution problems, 
mainly for second-order diffusion equations liisiiiaisa, but also for chemotaxis systems [7|, for 
non-local aggregation equations [nma, and for variants of the Boltzmann equation [20j . 

Lagrangian schemes for fourth order equations are relatively rare. Alternative Lagrangian 
discretizations for 0 or related thin-film type equations have been proposed and analyzed in 
[iiiini, but no rigorous convergence analysis has been carried out. We also mention two schemes 
iniiis] for the quantum drift diffusion equation, that is formally similar to (jlj) . In [53], the idea to 
enforce dissipation of two Lyapunov functionals has been developed, and was used to rigorously 
study the discrete-to-continuous limit. For the analysis here, we shall borrow various ideas from 
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A comment is in place on related Lagrangian schemes in higher spatial dimensions. Here, 
and also in our related works [STJ [28l [29] , the most significant benefit from working on a one¬ 
dimensional interval, is that the space of densities is flat with respect to the L^-Wasserstein 
metric; it is of non-positive curvature in higher dimensions, which makes the numerical approx¬ 
imation of the Wasserstein distance significantly more difficult. Just recently, a very promising 
approach for a truely structure-preserving discretization in higher space dimensions has been 
made [Sj. There, a numerical solver for second order drift diffusion equations with aggregation 
in multiple space dimensions is introduced that preserves — in addition to the Lagrangian and 
the gradient flow aspects — also “some geometry” of the optimal transport. These manifold 
structural properties enable the authors to rigorously perform a (partial) convergence analysis. 
It is currently unclear if that approach can be pushed further to deal with fourth order equations 
as well. 

Among the numerous non-Lagrangian approaches to numerical discretization of thin film 
equations, we are aware of two contributions |28L 134) in which the idea to enforce simultaneous 
dissipation of energy and entropy has been implemented. The discretization is performed using 
finite elements [53] and finite differences |53|, respectively. The discrete to continuous limit has 
been rigorously analyzed for both schemes. In difference to the convergence result presented 
here, certain positivity hypotheses on the limit solution are either assumed a priori |34j . or are 
incorporated in the weak form of the limit equation [23j . See, however, |2Ij for an improvement 
of the convergence result. 

The primary challenge in our convergence analysis is to carry out all estimates under no 
additional assumptions on the regularity of the limit solution Ui^. In particular, we would like 
to deal with compactly supported solutions of a priori low regularity at the edge of the support. 
Also, we allow very general initial conditions . Without sufhcient a priori smoothness, we 
cannot simply use Taylor approximations and the like to estimate the difference between {ua},. 
and Uif. Instead, we are forced to derive new a priori estimates directly from the scheme, using 
our two Lyapunov functionals. 

On the technical level, the main difficulty is that our scheme is fully discrete^ which means 
that we are working with spatial difference quotients instead of derivatives. Lacking a discrete 
chain rule, the derivation of the relevant estimates turns out to be much harder than for the 
original problem 0. For instance, we are able to prove a compactness estimate for maj but not 
for its inverse distribution function, although both estimates would be equivalent in a smooth 
setting. This forces us to switch back and forth between the original Q and the Lagrangian (101 
formulation of the thin-film equation. 

1.5. Key estimates. We give a very formal outline for the derivation of the two main a priori 
estimate on the fully discrete solutions. 

The first main estimate is related to the gradient flow structure of Q : it is the potential flow of 
the modified Dirichlet energy with respect to the Wasserstein metric >V 2 . The consequences, 
which are immediate from the abstract theory of gradient flows are that t {u(t)) is 

monotone, and that each solution “curve” 1 1 -^ u{t) is globally Holder-^-continuous with respect 
to W 2 . In order to inherit these properties to our discretization, the latter is constructed as a 
gradient flow of a flow potential (which approximates in a certain sense) with respect to 
a particular metric on the space of monotone vectors (which is related to W 2 ). See Section 2.1 
below for details. The corresponding fully discrete energy estimates are collected in Proposition 


11 


We are not able to give a meaning to the full energy dissipation relation (|^ on the discrete 
side, but this is irrelevant to our analysis. 

The second, equally important discrete estimate mimicks 0. Unfortunately, the L^-norm of 
Uxx is an inconvenient quantity to deal with, for two reasons. First, we need to perform most of 
the estimates in the Lagrangian picture, where 
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is algebraically more difficult to handle than the equivalent functional 

n p]\4 

/ {\ogu)l^dx = 


IQ 


IQ 




(19) 


which we shall eventually work with, see Lemma 12 Our discretization (12) is taylor-made in 


such a way that entropy dissipation yields a discrete version of (19). 

Second, the formulation of an iJ^-estimate would require a global C^d-interpolation of the 
piecewise constant densities ua that respects positivity, which seems impractical. Instead, we 
settle for a control on the total variation of the first derivative d^UA of a simple C'^’^-interpolation 
uAj see Proposition [T^ This TV-control is a perfect replacement for the iL^-estimate in ([^, and 
is the source for compactness, see Proposition [T7| 


1.6. Structure of the paper. Below, we start with a detailed description of our numerical 
scheme as a discrete Wasserstein-like gradient flow and discuss structural consistency of our 
approach. In Section we derive various a priori estimates on the fully discrete solutions. This 
leads to the main convergence result in Proposition |17[ showing the existence of a limit function 
It* for A —> 0. This limit function satisfies the weak formulation of ([^ stated in 0 ; this is 
shown in Section]^ Finally, we report on numerical experiments and discuss the observed rate 
of convergence in Section 


2. Definition of the fully discrete scheme 


The main aim of this section is to interprete the discrete equations (12) as time steps in the 
minimizing movement scheme for a suitable discretization of the functional with respect 
to a Wasserstein-like metric on a finite-dimensional submanifold Vsi^i) of 7^2 (f^)- 


2.1. Ansatz space and discrete entropy/information functionals. Fix K e N, let 6 := 

l/A, and define = Mk/K for /c = 0,1,..., A. For further reference, we introduce the sets of 
integer and half-integer indices 

= {0,1,...,A}, !+={!,...,A-1}, and = {1, , if _ ^ |. 

First and second order central difference operators DJ and are defined in the usual way: if 

y = {yi)t€iK ^ discrete function defined for integer indices f e 1/^ (i.e., on the nodes £^i), then 
DJ y and y are defined on half-integer indices n e (i.e., on the intervals Kk-I) Cre-i-i])) and 
on the “inner” integer indices fc e respectively, with 

rr-il-^1 rT'i2_-^n Uk+I ~ ‘^Vk + I/k-I 

yJk = g . yjfc = ^2 

If y = defined for half-integer indices A e instead, then these definitions are 

modified in the obvious way to have DJ y and y defined for integers k e Ik and half-integers 
K e , respectively; y needs to be augmented with additional values for y_i and in this 

case. 

Next, we introduce the set of monotone vectors 

y<5 := {(a;o, ■ • ■,xk) I To < ail < ... < xk-i < xk} ^ 

Each X e yi corresponds to a vector z = {zi/ 2 , Z 3 / 2 } ■ ■ ■ t ZK- 1 / 2 ) of density values Zk via 0- 
Our convention is that z_i = zi and z^^i = z^-i- For a function / : —> K, its first and 

second differential, d^f : y ^ and dif : y ^ respectively, are defined by 

[dsif{^)]k = dx^f(x) and by [dsf{x)]ki = d:^^dooJ{x). Further, /’s gradient V^/ is given by 
V 5 /(x) = 6~^d^f{x). For vectors v, w e the scalar product (•, •>^ is defined by 

^ /- 

<v, w >5 = 5 '^VkWk^ with induced norm ||v||^ = y (v, v)^. 

k=0 
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Example 4. Each component of z = zs[i] is a function on and 

d^z^, = -zi - 


( 20 ) 


where e M*" ^ is the kth canonical unit vector, with the convention eg = = 0. 

The main object of interest is the hnite-dimensional submanifold Vsi^) of 7^2 (f^) that consists 
of all locally constant density functions of the form u = U 5 [x], with given in (13), where 
X e To each density function u = U 5 [x] e Vs(S^) we associate its Lagrangian map as the 
monotonically increasing function X = X 5 [x] : [0,M] —> that is piecewise linear with respect 
to and satishes X(^fc) = Xk for k = 0,K. The density u and its Lagrangian 

map X are related by 

uoX = —. 

Remark 5. In one space dimension, the Wasserstein metric onV 2 {^) is isometrically equivalent 
to the L^-norm on the flat space of Lagrangian maps, see e.g. [33]. Our norm ||x —y||^ is 
not identical but equivalent to the L^-norm between the Lagrangian maps X 5 [x] and X 5 [y]. 
Consequently, there exist K-independent constants ci,C 2 > 0, such that 


Cl l|x-ylli < >V2(u5[x],U5[y]) ^ C 2 ||x-y||^ , for all x,yey5. 
See |5S1 Lemma 7] for a proof. 

Next, consider two functionals : y^ ^ K given as follows: 


( 21 ) 


H^^^(x) = hW(u5[x]) = f U5[x]log(u5[x]) dx = d 2 log(zfe_i), 

k=i 

h[^^(x) = H[2](u5[x]) = ^ / (u5[x])^da: = ^ XI 

Here is just the restriction of the the entropy to y^, and is the restriction of the 
quadratic Renyi entropy 

Zd«. 

Using (20), we obtain an expicit representation of the gradients. 






kgI 


1/2 


/ter 


1/2 


( 22 ) 


and — for further reference — also of the Hessians, 

d|HW(x) = d X 


kgI' 


1/2 


d|HL2J(x) = - X 


Kiel 


1/2 


(23) 


A key property of our simple discretization ansatz is the preservation of convexity. 
Lemma 6. The functionals and are convex on y. 


Proof. This follows by inspection of the Hessians (23). 


□ 


A conceptually different discretization is needed for the energy functional £ from ([^, which 
is identically +oo on Vs{Ll). 




HU- 
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Substitution of the explicit representations 
Es{x) = (WsUf\x),Wslif\x) 


in the definition (24) yields 

S 

~ 0 2-1 


kel 


1/2 


It remains to define a discrete counterpart for the potential V. A change of variables yields in 
the definition in © yields 

V(u) = / V{x)u{x) dx= V{X) dC, 

in Jo 

Thus, a natural discretization of V is given by 

V5(x) = d 2 V{xk). 

fceljf 

In summary, our discretization of is 

(x) = E 5 (x) + Vsix) = (WsU^^\x),Wslif\x)) + Vs{x). 


2.2. Discretization in time. Next, the spatially discrete gradient flow equation 

f = -ViE^ (x) (25) 

is discretized also in time, using minimizing movements. To this end, fix a time step with r > 0; 
we combine the spatial and temporal mesh widths in a single discretization parameter A = (r; (5). 
For each ye introduce the Yosida-regularized energy E^(-;y) : ^ K by 

EA(x;y) = ^ ||x-y||^ + Ey(x). 


A fully discrete approximation xa = (x^,x^,... ,x^,...) of (25) is now defined inductively from 
a given initial datum x^ by choosing each x^ as a global minimizer of E^(-; Below, we 

prove that such a minimizer always exists, see Lemma 

In practice, one wishes to define x^ as — preferably unique — solution of the Euler-Lagrange 
equations associated to EX(-;x^~^), which leads to the implicit Euler time stepping: 


-?n-l 

^ A 


= = - ^ (4hW(x) . d5^HM(x) + d|H["](x) . <5 jHW(x 


(26) 


Using (22) and (23), a straight-forward calculation shows that (26) is the precisely the numerical 
scheme (12) from the introduction. Equivalence of (26) and the minimization problem for E^ is 
guaranteed at least for sufficiently small t > 0. 

Proposition 7. For each discretization A and every initial condition x® e y^, the sequence of 
equations (26) can be solved inductively. Moreover, if t > 0 is sufficiently small with respect to 
6 and E^(x'^, then each equation (26) possesses a unique solution with E]^(x) sg E 5 (x°), and 
that solution is the unique global minimizer o/E^(-;5c^~^). 

Remark 8. In principle, the proof of Lemma^ below provides a criterion on the smallness of 
r > 0 that would guarantee the unique solvability of (26). We shall not make this criterion 
explicit, since in practice, we observe that the Newton method applied to (26) and initialized with 
always converges to “the right” solution x^, even for comparatively large steps t and in 
rather degenerate situations; we refer the reader to our numerical results in Section^ 

The proof of this proposition is a consequence of the following rather technical lemma. 

Lemma 9. Fix a spatial discretization parameter 6, and let C := £'^(ic°). Then for every y e y^ 
with E^ (y) ^ C, the following are true: 

• For each t > 0, the function E^(-;y) possesses at least one global minimizer x* G y^, 
and that x* satisfies the Euler-Lagrange equation 


= -WsEjm. 
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• There exists arc > 0 independent off sueh that for each r e (0, tc), the global minimizer 
X* G is strict and unique, and it is the only critical point o/E^(-;y) with £^(5) ^ C. 


Proof. Fix y e with (y) ^ C, and define the nonempty (since it contains y) sublevel set 

Ac := (E^(-,y)) ^([0,(7]) c y^. Let z = Z 5 [x], and observe that ^ 5/{b — a) for each z G 
From here, it follows further that 


Zk 


fc€l+ 





6 

Zi, I 1 + Zi,_i 
2 ^2 


1 

2 


U ^ (Zfc+1 + Zfc_i) 

V 



1 

2 




This implies that the differences — ^k.-^ = ^jz^ have a uniform positive lower bound on 

Ac- It follows that Ac is a compact subset in the interior of y^. Consequently, the continuous 
function E^(-;y) attains a global minimum at x* G y, 5 . Since x* G Ac lies in the interior of x, it 
satisfies difE^(x*;y) = 0, which is the Euler-Lagrange equation. This proves the first claim. 

Since E^ : y^ —> M is smooth, its restriction to the compact set Ac is A^-convex with some 
Ac e K, i.e., d|Ey (x) ^ AcIrt-i for all x g Ac- Independently of y, we have that 

d|EX(2,y) = diE,^(5) + ^llA-i, 

which means that x i—> EX(x, y) is strictly convex on Ac if 

5 

0 < CT < rc := , ■ ' ■ 

(-Ac) 

Consequently, each such EX(-, y) has at most one critical point x* in the interior of Ac, and this 
X* is necessarily a strict global minimizer. □ 


2.3. Spatial interpolations. Consider a fully discrete solution xa = (xai^a,---)- For no- 
tational simplification, we write the entries of the vectors Ica and ^ = Z 5 [x(^] as Xk and z^, 
respectively, whenever there is no ambiguity in the choice of A and the time step n. 

Recall that = u, 5 [x(^] g Vsi,^) defines a sequence of densitites on 12 which are piecewise 
constant with respect to the (non-uniform) grid {a,xi,... ,XK-i,b). To facilitate the study of 
convergence of weak derivatives, we introduce also piecewise affine interpolations zT: : [0, Ml ^ 
M>o and : 12 M>o. 

1 /2 

In addition to = kS for k G 1;^, introduce the intermediate points for k G . 

Accordingly, introduce the intermediate values for the vectors Xa and z^: 

i for k G 

Zfc = i(zfc+i -t Zfc_i) forfcGl+. 

Now define 

• zX : [0, M] —> M as the piecewise affine interpolation of the values (zi, Z3,..., Z;^_ i) 
with respect to the equidistant grid (^ 1 ,^ 2 ,...,and 

• Ma : 12 —> M as the piecewise affine function with 

uloXl = ^. (27) 

Our convention is that z^{f) = zi for 0 ^ ^ ^ 6/2 and Za(0 = ^ ~ '^/2 < ^ 

and accordingly u\{x) = zi for x e [a, a;i] and u\{x) = 1 for x G \xi^_i, b]. The definitions 

have been made such that 

Xfe = XA(Cfc), Zk = = u{xk) for all/c G Ik u 


(28) 
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Notice that 2^ is piecewise affine with respect to the “double grid” {xo,xi,xi,... ,Xj^_i,Xk), 
but in general not with respect to the subgrid (tq) , xk)- By direct calculation, we obtain 

for each fc e Ij- that 


d.ul 


d.ul 


(x i,Xk) 


(Xk,X l) 


Zk - 2fc_l 

Xk~Xk-k 
Zk+k - Zk 
Xk+k ~ Xk 


^fc+i 


Xk Xk—1 

Zj, I 1 — Zi,_i 

^fc+1 


= Zi._i 


fc+i ^k—h 


^k+Y 


Zu,l — Zu_l 

r\j ey r\j ey 


(29) 


Trivially, we also have that dxU vanishes identically on the intervals (a,xi) and (x^_i,&). 


2.4. A discrete Sobolev-type estimate. The following inequality plays a key role in our 
analysis. Recall the conventions that z_i = zi, Zj^_^_i = Zk-i, and that Zk = ^{zj.+ i + Zf.-i)- 


Lemma 10. For any x e is, 


5 Zk 
fcelt 


Ziui 1 — Zu 1 

A-l- o ^ o 


< - 

4 






1/2 


3 ( 4- 2^k,—1 

■ p 


(30) 


Proof. Define the left-hand side in p0| as (A). Then: 
(A) = 2-3 ^ Zfe(Zfe+l - Zfe_i )3(Zj,_^l - Z^_l) 


fcelj 

K 


S [zfc-l(Zfc_i - Zk-sf - Zk{Zk+i - Zk-0^] 

k = l 

S-^ ^ r 1 

Y ^k-k [Pk-k + )(^fe-i - Zk-kf - (Zfe+1 -f Zfe_l)(Zfc+l - 2fc_i)^J 


fc=l 


A 

> , Zi._i (zi._3 — Zi._i)(zi._i — Zi._3)3 -I- (Zk-l — Zk+l)(Zk + k — Zk-lY 


k = l 


+ 2zi,_ 1 (zi._i — Zk-i)^ — 2zi._ 1 (zi., 1 — Zi._i)3 . 

Rearranging terms yields 

(A) = -(A) - 2-3 Y 4-k [Pk-k - Zfc-3)3 - (zfe+i - ^fc-i)^] 


fc=i 


and further using the identity (a3 — P) = {a — b){af + P + ab), 
2-3 ^ 


G) --Vs 4-j [(= 


z^k-k -z, 




2^_|_ 1 — Zi^_ 1 


/c=l 


)“] 



k=l 


[D^Vfe-i[( 


Zu_i — -2l_ 3) + (z^ii — '2^_i) 

A. 2 ^ 2 ^ V ri^ 2 ^ 2 ^ 

-f (Zfc_i - Zfc_3)(Zfe+l - Zfe_l)j. 












CONVERGENCE OF A SCHEME FOR A THIN FILM EQUATION 

Invoke Holder’s inequality and the elementary estimate ab ^ + 6^) to conclude that 






k— i 




K 


s 4[DMd <-’2; 


Zk-k 


kgI 


1/2 


k=l 


(Zfe_i - + (; 


^fe+i Z.l,_l 


k— i 


)‘] 


( 31 ) 




where we have used an index shift and the conventions z_i = zi, Zj^j^i = z^_i in the last 
step □ 

3. A PRIORI ESTIMATES AND COMPACTNESS 

3.1. Energy and entropy dissipation. Fix some discretization parameters A = (r;(5). Be¬ 
low, we derive a priori bounds on fully discrete solutions (x^)®^g that are independent of A. 
Specifically, we shall prove two essential estimates: the first one is monotonicity of the energy 
E^, the second one is obtained from the dissipation of the auxiliary Lyapunov functional 
We begin with the classical energy estimate. 

Proposition 11. One has that Ej is monotone, i.e., Ej (x^) ^ E\{^ ^), and further: 

E^(xX) < (A) for all n > 0, (32) 

||x^ — x^ll^ ^ 2E^(x^) (n — n)r for all h ^ n > 0, (33) 


E 


^A 


= r2 ||V5E,^(xl)g<2E,^(4). 


(34) 


Proof. The monotonicity ( |36| ) follows (by induction on n) from the definition of x^ as minimizer 

ofEX(-;5^-i): 

^ si ^ A ) ^ ^||5X-^^-'0 + Er(5X) = EK(51;51-1) ^:EK(x1-1;x^- 1) =Er(x^-i). (35) 

Moreover, summation of these inequalities from n = n+ lton = n yields 


E 

n=n+l 






For n = 0 and n 00 , we obtain the first part of (34). The second part follows by (26). If 


instead we combine the estimate with Jensen’s inequality, we obtain 


1-^ 


II _*n II , V”' 

n=n+l 


^ r 


E 


n=n+l *- 


||.^ 1|| 


2 \ 1/2 

j (T(h-n))^''/ 


which leads to (33). 


□ 


The previous estimates were completely general. The following estimate is very particular for 
the problem at hand. 

Lemma 12. One has that Hjj^^ is monotone, i.e., H W(xX) < hW( x^ ^). Moreover, it holds 
for any T > 0 that 


N 


- E ^ E ("")' 

"=0 «e4" 
for each Nj. e N with Nt.t e (T, T + 1). 


3 / + z^_i 


< 


4(H[i] + AM(T-H)), 


<52 


(36) 
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Proof. Convexity of implies that 


hW(5^-1)_hW(x1) > (v5HW(xX), 21-' = r(v^HW(5X),V^Er(5X)/^, 


for each n = 1,..., N^.. Summation of these inequalities over n yield 


N.r 


Y , (v^hW(^^), V5E,^(x1)) ^ hW(4) - hW(x^). 


(37) 

To estimate the right-hand side in (37), observe that Hj^^^(x^) ^ Hbl by hypothesis, and that 
(x^) is bounded from below thanks to the convexity of s sln(s) and Jensen’s inequality, 
which yieds for any x G 

Hj^^^(x) = J U 5 [x] lnu 5 [x] dx > Min ■ 

We turn to estimate the left-hand side in (37) from below. Recall that E^ = E^ -f V^. For 
the component corresponding to V^, we find, using (22) and 

'v,H["(5ri),v,v,(xi)> -i 




1/2 


:> 


infF,;,j(x)) d V 
£ceR I ^ 


> -AM. 


ftel 


1/2 


The component corresponding to E^ is more difficult to estimate. Thanks to (22)&(23), we have 
that 

[1], 


4^V5E5(x),V5HL^J(x 
= 4 /V^H^ (2), (2) V^hW (x)\ + 4 /v^H^ (2), (2) (x)' 


= 2J 2 


K€l 





K€l 



Further estimates are needed to control the second sum from below. Observing that 




— 2z'f + z'^_i — 2zk + Zi^_i 

- = ZZk' --2^-- + 




d2 


Z^+l - 


Z„_1 - 


and that 2ab > ~ for arbitrary real numbers a, b, we conclude that 


4^V5E5(x),ViHW(2)^^ 


> 4 






1/2 


+ Zi^—i 


2(5 

T 2. 


Keljc 


-Z 


fc-i 


Now apply inequality (30). 


□ 


3.2. Compactness. The following lemma contains the key estimate to derive compactness of 
fully discrete solutions in the limit A —> 0. Below, we prove that from the entropy dissipation 
(36), we obtain a control on the total variation of dxu\. 

Several equivalent definitions of the total variation of / G L^(f2) exist. In case of piecewise 
smooth functions with jump discontinuities, the most appropriate definition is 

TV [/] = sup I ^ \f{rj+i) - f{rj)\ : J G N, a < ro < ra < • • • < rj < . (38) 
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Further recall the notation 


Ifh = lim/(a;) -lim/(a:). 

xlx XJX 


for the height of the jump in /(a;)’s value at x = x. 

Proposition 13. For any T > 0 and Nr e N with xNr e (T, T + 1), one has 


Nr 


r ^ TV ^ Ct. 


(39) 


Remark 14. The proof below yields for Ct the explicit value 


Ct = 50(H[i] + AM(T+ 1)). 

Proof. Fix n. The function dxuf^ is locally constant on each interval {xj._i, Xk), and equal to 
zero elsewhere. Therefore, the total variation of is given by the sum over all jumps at the 
points of discontinuity, 

TV[d,21]= 2 \ldxulU\+ E lld.ulhJ- (40) 

fc€l+ 

The jumps can be evaluated by direct calculation: 


Zk-Zf,_l Zfe+I-Zfe.l 


Xk-x,,_i Xk-Xk-i ^ ’ 


dxU 

This implies that 


_ Zk+i - Zk _ Zfc+i - Zfc_i _ 
(xk.x^^i) Xk+i-Xk Xk+l-Xk 


Zu, 1 — Zu_ 1 

• V n III n 


Zuii Zu _i 


for k G 


for k G 


(41) 


II<5. 


= <5 


fc+i k—-k 


for fc G Ij-, 


II<5: 


X ‘^A NIEk; 




j_ 

P 


for K G I 


1/2 

K ■ 


We substitute this into ( |40[ ), use Holder’s inequality, and apply (301 to obtain as a consequence 


of elementary estimates that 
TV [dxul\ ^ 6 Y, 


\ 2 

^fc+i ^k-i 


feel! 


1/2 


< 




Si 

vfcelt 


feelt 


KEl 


' lyU _ .yU 

/c+i /c—i 


4 \ 2 


+ u (4)"[D^-i]^ 


keT 


1/2 


;(2(6-a))i is Y 


K€l 


1/2 


We take both sides to the square, multiply by r, and sum over n = 0,..., Nr- An application of 
the entropy dissipation inequality (36) yields the desired bound (39). □ 
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3.3. Convergence of time interpolants. 

Lemma 15. There is a constant C > 0 just dependent on E'^ and (b — a), such that the following 
estimates hold uniformly as A ^ 0; The functions {ua}^. and {Sa}^ are uniformly bounded, and 

(42) 


Proof. For each n e N, 


sup 114 {SaIt (i)llL 2 (n) < C, 

sup II {ua}.^ (t) - {ma}^ {t)\\m(n) < CS, 

sup II {uA}r (i)l|L“°(n) < C. 


(43) 

(44) 


feel+ ^ fc+i fc 

_ ■yi ^K-i ~ / ^fc+l 

“ ^ 2 V a;” 


k— i 


™n _ ™n 


k-h\^ Xk-X^-l ( ^k+l 

2 V xl-xl_^ 




_i_ 

4+i + 4- 




4+i 


^2E5(xX). 


fcelt 


This gives (42|). For proving (|43|, we start with the elementary observation that 
Therefore, 


lulix) - UAix)\ ^ \z^^i -^4il for all xe [x^_i,x^+i]. 

1/2 


\\ul-ul\\LHn)^S Y, S 

keif 


~ » fc ~ » fc 

4+i 4-i 


1/2 


fc€l+ ^ \ keif 


4+i 4-i 


^<5(2(&-a))l/2E,(xX)l/^ 

which shows ( |4^ . Finally, ( [44| ) is a consequence of (42)fc([4^. First, note that 

||{2A}^(i)|lii(n) < II {uA}r WIIlTO) + II {uA}rit) “ {^A (i) || Li (O) ^ M + C6 

is uniformly bounded. Now apply the interpolation inequality 

II {Sa}^ WIU«>(n) < (^114 {Sa}^ WIIl2E)II WIIliE) 

to obtain the uniform bound in (44). 


□ 


Proposition 16. There exists a function : Kg,o x ^ K;>o that satisfies for any T > 0 

e C^/\[0,T]-,V2m n L”([0, T]; Ffi(f2)). (45) 

Furthermore, there exists a subsequence of A (still denoted by A), such that the following are 


true: 


{uaIt- (t) —> u,f{t) in V 2 {^), uniformly with respect to time, (46) 

{SaIt- —> uniformly on [0,r] x 12, (47) 

{XaIt- (t) —> X,|,(t) in L^([0,M]), uniformly with respect to t e [0,T], (48) 

where X,,, e C^/^([0, T]; L^([0,M])) is the Lagrangian map ofu,^. 


Proof. From the discrete energy inequality (33) and the equivalence (21) of W 2 with the usual 
L^-Wasserstein metric yV 2 j it follows by elementary considerations that 

v2 


W 2 ({uA}^(t),{wA}^(s)) ^C(|t-s|), 


(49) 


for all t,s e [0,T]. Hence the generalized version of the Arzela-Ascoli theorem from [U Propo¬ 
sition 3.3.1] is applicable and yields the convergence of a subsequence of ({uaIt) to a limit 
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in V 2 {^), locally uniformly with respect to t e [0, oo). The Holder-type estimate (49) implies 
u e C'^/^([0 ,oo); 7^2(^^))- The claim (48) is a consequence of the equivalence between the Wasser- 
stein metric on 7^2(^) and the L^-metric on X, see Remark]^ In addition, the limit function 
is bounded on [0,T] x H, thanks to (44). 

As an intermediate step towards proving uniform convergence of {ua},-) show that 

WA(i) —>■ u^it) in uniformly in < e [0,T]. (50) 

For t e [0,T], we expand the L^-norm as follows: 


{uA}r{t) [({wA}^-M=f) {uA}^](t,a:)da; 

' {uA}r - W*) ( {Sa}^ - {uA}r )] x) dx 
{uA}r - U=I=)m=i= j (t, x) dx. 


On the one hand, observe that 


sup 
te[0.T] 7 a 


{ma}^-«=!=)( {wa}.^ - {ua}^) {t,x)dx 


< sup ((II {ua}^ (i)IU®(a) + l|M*(t)||L<»(a))ll {ma}^ {t) - {ua}^ WIIlho)) 

ie[0.T] 

which converges to zero as A ^ 0, using both conclusions from Lemma On the other hand, 
we can use a change of variables to write 


[uA}r - ■“*) {ma}^ j {t, x) dx • 


UAlr ~ 

M 


)M,i,j(t,x) 


dx 


= J [{uA}T--u=i=](i,{XA},-(t,x))d^-y |^{2 a}t-- u=i=](i,X,f(t,5))d^. 

We regroup terms under the integrals and use the triangle inequality. For the first term, we 
obtain 

f-M 

({2a}. (i,{XA}. {t,0) - {2a}. de 


sup 

te [O.T] 


r-M /.{XA}^(t,5) 


^ sup / / 

te[O.T]7o 


< 


< 


\dx{uA}r \ {t,y) dyd^ 

pivi 

sup / II {ma}. ||Hi(a)|X,f - {Xa}. 
e[o.t] 7o 

sup (^11 {ua}. (<)||i/i(a)||X*(f) - {Xa}. (t)|| 

Z,2([o,M])) • 


te[0.T] 

A similar reasoning applies to the integral involving in place of {ma}.- Together, this proves 
(50), and it further proves that e L® ([0, T]; iL^(H)), since the uniform bound on 2 a from 
(42) is inherited by the limit. 


Now the Gagliardo-Nirenberg inequality (88) provides the estimate 

II2/3 


II {2a}. (t) - w=f(t)||ci/6(a) < CII {2 a}. {t) - MH=(i)||^((n)|| {2 a}. {t) - (51) 

Combining the convergence in by ( [50| with the boundedness in H^{n) from (42), it 

readily follows that 2a (i) ^ UH=(t) in C'^/®(H), uniformly in t e [0,T]. This clearly implies that 
{2a}. —> uniformly on [0,T] x H. □ 

Proposition 17. In the setting of Proposition [7^ we have that e L^^(R^o', o.'nd 

{2a}. ^ strongly in L^{[d,T\,H^{^l)) (52) 

for any T > 0 as A ^ 0. 
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Proof. Fix [0,r] c R^q. Remember that is differentiable with local constant derivatives 
on any interval for k e Ij- u u {K}, and it especially holds dxU^{x) = 0 for all 

X G (a, a + 5/2) and all x G (6 — 5/2, b). Therefore, integration by parts and a rearrangement of 
the terms yields 


ll^xMAlli2(n) = Yi j dxv5ldxu\dx= ^ ^u\{x)dxu\{,x)^^_^ 

< l|21llz,co(n) TV[d,21]. 


Take further two arbitrary discretizations Ai, A 2 and apply the above result on the difference 
{Sai},- — {wAalr- Using that TV [/ — ff] TV [/] + TV [g] we obtain by integration w.r.t. time 
that 


j-T ^ 

J \\5x {uAi}x ~ {UA 2 }r\\L^(Q) 


sSrU2 sup \\{uaJx- 

ie[0,T] 




2^ TV[d,{2Ajj" + TV[d,{SAaU" 


1/2 


This shows that {ua}^ is a Cauchy-sequence in L‘^{[0,T]-H^{n)) —remember (39) and especially 
the convergence result in (471 — and its limit has to coincide with u,,, in the sense of distributions, 
due to the uniform convergence of {ua}^ to Mh= on [0,T] x ft. □ 


4. Weak formulation of the limit equation 

In the continuous theory a suitable weak formulation for 0 is attained by applying purely 
variational methods, see for instance nzinn]. More precisely, the weak formulation in 0 is 
obtained by studying the variation of the entropy S along a Wasserstein gradient flow generated 
by an arbitrary spatial test function p, which describes a transport along the velocity field 
p'. The corresponding entropy functional is $(u) = /jj p(x)u(x) dx. It is therefore obvious to 
adapt this idea - similar as in [28l [29j - to show that {ua}t- inherits a discrete analogue to the 
weak formulation ( [T7| ). Hence, we study the variations of the entropy along the vector field 
generated by the potential 

$(x)=/ p(X^[x])de 
Jo 

for any arbitrary smooth test function p G with p'{a) = p'{h) = 0. That is why we define 

rM 

v(x) = V54 >(x), where [d 5 f>(x)]^ = / p'(X(^))6»fc(^) d^, k = l,...,K-l. (53) 

Jo 

Later on, we will use the compactness results from section [T3| to pass to the limit, which yields 
the weak formulation of our main result in Theorem Therefore, the aim of this section is to 
show the following: 


Proposition 18. For every p G C'®(H) with p'{a) = p'{h) = 0, and for every p G C'®(R>o), the 
limit curve satisfies 


rCC p pec 

/ / dtipu:fdtdx+ / N{u:f,ip)dt = 0, 

Jo Jn Jo 


where the highly nonlinear term N from (18) is given by 


N{u,p)=^ [ {u'^)xp'" + 3 uIp" dx + [ Vxup'dx. 
2 Jn Jn 


(54) 


(55) 


The proof of this statement will be treated in two essential steps 

(1) Show the validity of a discrete weak formulation for {ua},., using a discrete flow inter¬ 
change estimate 

(2) Passing to the limit using Proposition 


17 
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For definiteness, fix a spatial test function p e with p'{a) = p'{h) = 0, and a temporal 

test function rj e C®(K>o) with supp?? c (0,T) for a suitable T > 0. Denote again hy Nr e N 
an integer with rNr e (T, T + 1). Let B > 0 he chosen such that 

WpWcHQ) < B and hllcbRso) < B. (56) 

For convenience, we assume i5 < 1 and r < 1. In the estimates that follow, the non-explicity 
constants possibly depend on T, B, and E^, but not on A. 

Lemma 19 (discrete weak formulation). For any functions p e (7®(n) with p'{a) = p'{h) = 0, 
and rj G C'”(K>o), the solution x)^ with u\ = U 5 [x^] of the minimization problem ( |15[ ) fulfills 
00 

n=0 

where we use the short-hand notation p'{x) := (p'(a;i),... ,p'{xk-i)) for any xG 

Proof. As a first step, we prove that both vectors p'{x) and v(x) nearby coinside for any x G y^, 
i.e. it holds 


^(x^) ^^(x^ ) _ (xl), = 0{r) + (57) 


||v(x)-p'(x )||_5 (58) 

for a constant C > 0 that only depends on B and ft. Hence, denote by X = X 5 [x] the corre¬ 
sponding Lagrangian map of x and choose any fc = 0,..., iL, then one easily gets 

|v(xfc) - p{xk)\ = |[V5^(x)]^ - p{xk)\ ^ / |p (X(^))6'fc(^) - p 

Jq 

First assume ^ G [^fe_i,^fc], then a Taylor expansion for p'(X(^)) yields 

p'(X(e)) = p\xk) - - a) = p'{xk) - /(X(a))(aifc - a:fc_i)0fc_i(C) 

for a certain G Consequently the valitdity of ^ OuOk-i d^ = ^ yields 

|p'(X(a)0fc(C) - p'{xk)9k{0\di ^ lB{xk - Xk-i). 
hk-i ^ 

Similarly one proves the analogue statement for ^ G hence 

3 3 

|v(a;fc) - p'{xk)\ < ^{{xk - Xk-i) + {xk+i - Xk)) = -B{xk+i - Xk-i). 

Squaring the above term and summing-up over all fc = 1,..., AT — 1 finally proves (58), due to 
{xk+i - Xk-i) < 2(6 - a) and 

qa ^—1 qc ic—i 

||v(x) - p'{x)\\l < —B"^ X i^k+i - Xk-i)"^ -TxB'^ib - a) X (^k+i - Xk-i) < C5. 

k=l k=l 

Let us now invoke the proof of ( [F^ . a Taylor expansion of p for X,X' G X yields 

p(X) - p(X) - |(X - X')^ < p'(X)(X - X), 
which implies for X' = X 5 [x)^“^] and X = X 5 [x)^] 


$(xl)-$(xl-i) 


B 


-fC-l nM 

p'(x,[xi])0fc(c)de 

fe=i Jo 

= <x^ - x^-1, V,#(xl)X = r (5^), V,$(xl)X • 

(59) 


Thanks to ( [5^ , the last term can be estimated as follows: 

r <V^E,^(^1), Vi$(xX)>, = r <V^E,^(5^), p'(^)X + ^ (V^E^(^), v(xX) - p'(51)>, 

< r<V5E,^(xl),p'(x^)X + r ||V5E,^(xl)||, ||v(x^) - p'(xl)||, (60) 
^ r <ViE,^(21),P (51)X + Ct5^I^ ||ViE,^(xl)||, . 
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So combine ( [59| and (59), add r]{{n — 1)t) and summing-up over n = I,^ N, one attains 


r Y, V{{n-1)T) 


71=1 




-<V^E,^(xX),/9'(S1)X 


<II^IIc’ 0 ([o.t])T XI ^ 11^ “ + C'll^llc<)([o.r])T X (^)\\s ’ 

n=l 71=1 

where the right hand side is of order 0{t) -f 0(5^/^), due to ( [M| ). An analog calculation replacing 
p with —p leads finally to (57). □ 


The identihcation of the weak formulation in (54) with the limit of (57) is splitted in two main 
steps: In the first one, we estimate the term that more or less describes the error that is caused 
by approximating the time derivative in (54) with the respective difference quotient in (57), 

rT 


SiA 


/o 


m 


V'{t) l^p{x){uA}r{t,x)dx + p{t){(Ws'El{xl),p'{5i^))g}^{t)j dt 


(61) 




The second much more challenging step is to prove the error estimate 


62 ,A 


vit) 


1 


p"'{x)doo{{uA}l){t,x) + ^p"{x)d^ {ua} 1 {t,x) da; 


+ J^V^{x){u a} rP{x) da; - {<V5E^(5c2,), P (xl)> (t) j dt 


(62) 




which, heuristically spoken, gives a rate of convergence of {(V^E^(x^), p'(5c^))>^}^ towards the 
nonlinear term N(u^,p) from (55). 

The Hrst estimate in (61) is a consequence of Lemma [I^ 

Proof of (61). Using that pinr) = 0 for any n ^ N^., we obtain after “summation by parts”: 

- [ v'{t) ( [ p(a:){uA}.x(^>a^)da;) dt = - X ( [ v'{t)dt [ p(a;)MA(a;) da; J 
Jo Van / \a(n-i)T an J 


Ax 


--X 

71=1 \ 

( 


( pjriT) - 7?((n - 1 )t) 


rM 


poXliOdf 


(63) 


Nr / pM 

X X 1 )'^) / 

n=l V -^0 


pox^(a-poxr'(0 


de 


Finally observe that 

rT 


R := 


Ax 


v{t) {<p'(xa), V^E^(xa)X}^ (t) dt - T X v{{n - 1)t) (V^E^(x^), p'(x^)>_, 

^6 n=l 

1/2 


Ax 

s: I T X 

71 = 1 


[ p{t) dt - p{{n - 1)t) ] [7 

a(n-l)T J \ 


1/2 


X ||V5E,^(x1)| 


71=1 


s: ((T-f- 1)BV2)^/^(2B2 e^(j4))1/2 = C"Ef^{^ aY/'^t, 
using the energy estimate ( |34| . We conclude that 

Ax / 


ei.A ?i? + rX (M(n-l)r)| T _ <v^Er(^), p'(51)> 

1=1 \ Jo T 


= o{t) + 
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where we have used (57), keeping in mind that ^(x^) = pO^a) 


□ 


The proof of (62) is treated essentially in 2 steps. In the first one we rewrite the term 
(V5 EJ"(x 2^), /9'(5cjy)^ (see Lemma 20), and use Taylor expansions to identify it with the corre¬ 
sponding integral terms of (55) up to some additional error terms, see Lemmata 23]|^ Then we 
use the strong compactness result of Proposition to pass to the limit as A ^ 0 in the second 
step. 

Lemma 20. With the short-hand notation p'(x) = (p'{xi),..., p' {xk-i)) for any x e one 
has that 


- (Vs-Ej{^A),p'm}s = A-,+A^ + A^-A^ + A^,+A^ + A^, 


(64) 


where 

A- = 6Y, 


kelt 




feel; 


■^5 = 3 S 

feel); 

^4=^2 

feel+ 

^5=^2 


^fc+i ^fc-i 




^k+h ^fe-i 


rr-ik 

^fc+i ^k-i 


n 32 I 32 , n n 

+[Z,_i) +Z, iZ,_i 


p'(x^+i) - p' 


2S 






p'ixl^^) - p'{x^) 


p'{xl) - p’{xl_^) 
6 

n \3 \ 


feel 


^6=^2 


feel 



22 " iz" 1 

^ fe+2 fe-2 


p"i^k), 

p'{x^k+i)-p'K)-K+i-^k)p''K) 

(52 


(^fc+i)' + fp'{xl_,) - p'{xl) - {xl_^ - xl)p"{xl) 


(52 


A? = ^ 2 V{xl)p\xl). 


feei: 


Proof. Fix some time index n G N (omitted in the calculations below). Recall the representation 
of as 


1 


V^E^S) = ^ (d|HW(^)^?.HL^J(^) + 5|HL^J(x)djHL^J( 


rWi 


with corresponding gradients and hessians in (22) and (23). Multiplication with p'(xa) then 
yields 


■<V5Er(5^),p'(kA)>, = I 2 zl 


1/2 


-k+4 


1 -2z^ + Z^_i\ f p'{x^+i) - p'{x^_i)' 


K6l 




(52 


- 2^2 + z2_ \ / - p’{x^l)' 


«el' 


1/2 


(52 


+ ^ 2 V{xk)p'{xk). 


feei; 
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Observing that 


zj+i ~ + zi_^ 

52 

we further obtain that 


= 2z, 


z„+i - 2z„ + z„_i ( z^+i - Zk 


52 


+ 


Zk -1 - Zf, 


-(VsEy(x^),p'(x^)\ = S z^f 


z^+i - 2 z«; + z^_l^^ / p'(x^_^i) - p'(x^_i)' 

<52 


+ A .2 + -^43 + ^7. 

It hence remains to show that (A) = Ai — A 4 Aq^ where 

'Z^+I -2z^ + z^_i\ (p'{x^^i)- p'{x^_iy 


{A) 4 2 


kgI 


1/2 


<52 


After “summation by parts” and an application of the elementary equality (for arbitrary numbers 
p± and q±) 

P++P-, ^ , s 9+ + 9- 

P+q+-P-q- = -;;-( 9 + - 9-) + [P+ -P-) --, 


one attains 


o) = 5E 

fc€l+ 

2 2 

fcelt 


^k+A - Zk-1\ f^k-k ^k+k\ f p'{Xk+l) - p'ixk-l) 


^fc+i - ^fe-i \ f^k-k + ^k+k \ f p'{xk+i) - 2p'{xk) + p'ixk-i) 




kelt 


^k+k - Zk-k \ l^k-k + ^fc+i '\ f p'{xk+i) - 2p'{xk) + p'{xk-i) 


, (65) 


where we additionally used the identity {p^ — = (p — q)(p^ + 9^ +P9) in the last step. In order 

to see that the last sum in (651 equals to —A 4 + A 5 + Ag, simply observe that the identity 

1 1 

Zhj-l Zu_i 

2^2 _ _ 
makes the coefficient of p"{xk) vanish. □ 


^fc + 1 ^k — 1 ^k 

j 5 


— Zi._i 
.^TJL _^ , 2,. 

Zh._i_ 1 Zu 1 

fi-i- 5 fi, 5 


For the analysis of the terms in (64), we need some sophisticated estimates presented in 
the following two lemmata. The hrst one gives a control on the oscillation of the z-values at 
neighboring grid points: 


Lemma 21. For any p,q ^ {1, 2} with p F q ^ 3 one has that 

p 


Nr 

"=1 fcelt 


fe+i ^fe-i 






1 


< 




( 66 ) 


Proof. Instead of ( 661 , we are going to prove that 

^ E ^ E ^ 

22=1 keif 


fe+i 






(67) 


is satisfied for any p,q e {1,2} with p + q ^ 3, which implies (661 because of the following 
considerations: The situation is clear for q = 1, thus assume g = 2 in (66). Then (67) is an upper 
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bound on (66), due to 




1 = 


+2 ^-1 


.^k+i 




for any n = 1,..., Nr ■ 

0 p: 

Nr 


To prove (67), we first apply Holder’s inequality, 

p 




K 


^10 ^10 
^k+k ^k-l 


"=1 fcel+ 
Nr 

fcelt 


k±k _ ^ 

k+k 


Nr 


^k+k ^k-k 


V P + 9 
4\ 4 


”=1 feel+ 
Nr 


^k+k ^k-l 


p+q / ^ \ ^ 

I 0 \ 


,h+k 


( 68 ) 


"=lfc€l+ V fc+2/ 

with a = 1 — . The first factor is uniformly bounded due to (30) and (36). For the second 

term, we use (44) and (891 to achieve 

. V 9 

Nr / e \ A 


E 

fc€l+ 


'fc+i , 


^ (T + l)S{b — a) “ II {m}^ ||L“o([o,T]xa) ^ C!{T + 1)(5(6 — a) ■ 


which shows (671, due to a > |. 

Lemma 22. For any p e {1, 2} one obtains that 


Nr 

^ H < S 4‘ 

”=1 fcel+ 

Proof. Appling Holder’s inequality, 

Nr 

^ 2«2 4 


_ ^IL 

k+h k-k 




□ 


(69) 


n=l 




/ Ij ^ 

k+h k-k 


Nr 


E^E 


ss I T j>, 0 } . z',: 
"=1 feelt 


^11 ^11 
^k+k ^k-l 


- xl_,Y 

4XV2 / 


1/2 


^ E E ^k{x 

"=1 keJk 


k+1 


rV- 




The first sum is uniformly bounded thanks to (30) and (36), and the second one satisfies 

1/2 


Nr 


■ 2 ■* 2 4(a 

^=1 kelt 


k+1 


I < snr + l)^/^|| {2}. lli/^([0,T]xa)(& -«) 


where we used (44) and (89). 


□ 


Lemma 23. There is a constant Ci > 0 expressible in O, T, B and FN such that 


Rk-.= rY, 


n=l 


Af-3 ^AiOdi^mrp"o^mdf 

Jo 


^ C'i(5++ 


Proof. Let us introduce the term 

Sr := 2 2 4 

fcelt 


^k+k ^k-l 




p"ox^(0 de 


2 
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First observe that by definition of 2^, 

pM 


^fe+i ^fe-i 




^0 

hence we get for 

nM 


fcelt 


nlVl 

3/ ?Sm21(0VoXl(0de-SI 
io 


^3S 2 


^fc+i ^fc-J 




/celt 


fcelt 


^fc+i ^fc-i 


'fc+i 


- 1 


fe— i 


This especially implies, due to (661 that 

Nr 




oM 


B^-3 22(6d52X(aV'oXX(Ode 


/o 




(70) 


For simplification of i?i, let us fix n (omitted in the following), and introduce e \xkTXk+i\ 
and x~^ G \xk-iTXk\ such that 

p'{xk+i) - p'{xk-i) _ p'{xk+i) - p'{xi^ ^ t[{xk) - p'{xk-i) 


25 


Recalling that 


= P”{it) 


25 

+ \ ^fc+1 


25 


+ P'\it) 


25 

Xk+i-Xk 1 (p''{.xl) p"(ife) 


25 




z^_i 


r^k+1 


' ^fc-i 


0fe(e) de = ^, 


(71) 


one has for each fc G I J- , 




P''{K) , P"(5fe) 


^fc+i 


‘fe-i 


7^) 




Zkp" o Xa d^ 




2 + 


Zk-l 


Zi-.i 


+ 7^fe+l 2 + 


fc+2 _ 
^fc+i' 


p”{x^) + 7^/c+i 2 + 


fe— i 




fc+5 


Zk_i 


p''(^k ) + 2 + 


'fc+1 


'fc-1 




r4fe+i 


z/cp" o Xa d^. 




and furthermore 

(^)=7 


^k-h 


Zi._i 

K 2 


-1 + z, 


fc+i 


+ 


fc-J 

z2 
/ ^2 




"fc+i 


1 +z 


/c+4 


fc— -4 


k— 4 


P"iik) 




. ..... 

Zfc [p” O Xa - p"(5fc )] ^ [p" o Xa - p"(xfe )] d^. 


Applying the trivial identity (for arbitrary numbers p and q) 




P 


q 
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the above term finally reads as 


(Al) = 


Zk-i ( ^-1 I + 2zfc I ^-1 


Zk+2 






^/c+i 


h. r, 


1 + 2zh 


^fe+i 




(72) 


3 


J Zk[p" oXa- p"{x+)]dC-^ Zk[p" oXa- p"{x^^)]dC 




Since lies between the values Xk and and Xa(C) s \xk,Xf.j^i^ for each ^ e [^fe,i], we 

conclude that |Xa(C) ~ I ^ Xk+i — Xk, and therefore 

^ , „ . .. 3 

2<5. 




Zk\p oXa(^)-P {xl)\diik -Bzk{xk+i-Xk). 


(73) 


A similar estimate is valid for the other integral over [^fc_ i, ^fe] and for the integrals with p"{x ^.). 
Thus, combining (72) and (73) with < 2^^ and the definition of A", one attains that 

2 r / -.n \ / -.n n 


\A- - B-\ < 2B 2 z," 

fc€l+ 

+ 35 X 2 




^fc+i i 




-1 + 


“fc -1 


- 1 




^fc +1 


k+h 


n \ 
fc-lb 


kelt 


and further, applying ( 66 ) and (^69 




T ^ |A^ -5^*1 ^ 


By triangle inequality, (70) and (^74| provide the claim. 


(74) 

□ 


Along the same lines, one proves the analogous estimate for A 2 and A 3 in place of A": 
Lemma 24. There are constants C 2 > 0 and C 3 > 0 expressible in 17, T, B and such that 




R2-.= rY, 

n = ] 

53 := r 2 


1 


nM 




n = l 


nM 


A? - 






Lemma 25. There is a constant C 4 > 0 expressible in 17, T, B and E^ sttc/i that 


54 := r ^ 


AJ- 




^ (74(5^/^ 


Proof. The proof is almost identical to the one for Lemma above. As before, we introduce 
the term 


S 4 ” := <5 




2 j/ 


fc€l+ 




/ox^(e)de 


and get due to ( 66 ), analogously to (70), that 

N.r rM 


•E 

n=l 


51 


0 








(75) 
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By writing 


one obtains that 


rV r\ rL ry 


9 


r\j fy rv fy 




„ ^ P'\^k) - T / Zkp" O X^(C) 

- l) + - l) - J 4[P" oXl(C) - /(o;^)] de 


Observing — in analogy to (731 — that 


1 


rik+i 


" ziv oXKC) - p"(x^)| dc < - ^Li) 




we obtain the same bound on \A 2 — as before on |A" — i?"|, i.e. 

\ 2 


IA2 - B2I ^bY,4 


kelZ 


'fc+4 


k— i 


1 + 


fc— 4 


fe+i 


+ 




’ ~f£. ~/£. 

^fc+i ^k-i 


(a^fe+i -a;fe_i). 


feel; 


Again, applying ( 66 ) and (69), we get 


n=l 


and the estimates (75) and (76) imply the desired bound on R 4 . 


(76) 

□ 


Lemma 26. There is a constant Cs > 0 expressible in 17, T, B and E'^ such that 




^5 := r 2 


- 1 ^liOd^^AiOp'" o xlio d^ 


< 




Proof. The idea of the proof is the same as in the previous proofs. Let us define similar to i?" 
the term 

r«fe+i 


:= <5 2 


^k+ i ^k -i 


"fc+i 


feei; 


Ys. 


2 Jff 


^k-l 


p'"oxi{0df. 


Note in particular that we weight the integral here with z'^ i. Then 


rM 


zl{OSizl{Op'"oXliOdf-B^ 


< 


>1. 


fe€i; 


< 


fe +2 ■ 

( .yU _ .yU 

fc+5 k—k 

Ibs Z 

fcelt 


rik+l 


Bk-l 






'^k— k 


- 1 


'■/c+i 


where we used that by definition of 


nM 


^iOSiZ^iOp"'o^AiOdi= 2 

fcelt 


y,IO _ y,IO 

^fe+i ^fe-i 


'^k- 
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This especially implies, due to ( 66 ) that 


n = l 






/o 




(77) 


Furthermore, one can introduce intermediate values such that 

1 
2 


- xl)p"{xl) = hxl^, - xlfp"\xt) = ,^J^p'"ii+). 


Using the identity 


1 


^10 ^10 

^+9 _j_ ^ 2 




we 


V 2 J 2{z-^,r 2 

thus have — using again ( |7T| ) — that 
/ i4+^f + - p'(x^) - - z^)p"(ai^) 


1 + 




fc+i / fc—i 


fc+4 


- 1 


"fc-l 




- 1 + 


^fe+i / ^fe -1 




(52 


/ 2(5. 


\ 1 

\ _ _j_ I ^ jn ^ vn 

) 25 4_, 


z^/'oX^dC 




Observing — in analogy to (73) — that 

K, 


1 


/ ^71 _// ^ -\rn (c\ ff f ^n\ J r' ^ ^71 /_71 _n \ 

-^fc |P ° — p (sj,)! d^ ^ —Zj,_|_i (a;j,_|_i — x^_i), 


and ^ 2 z^, we obtain the following bound on jAg — Bg |: 


B 




~»t- ~»t- 

^fc+i ^fe -1 


feelt 


B S 

fcelt 


~/£. _ ~/£. 

^fe+i ^fc-i 






-1 + 


ri 

fc-lJ- 




- 1 


k+h 


Again, applying ( 66 ) and (691, we get 




(78) 


□ 


n=l 

and the estimates and ( [T^ imply the desired bound on B 5 . 

Arguing like in the previous proof, one shows the analogous estimate for Ag in place of Ag . It 
remains to analyze the potential term A”, where we instantaneously identify the ^-integral with 
the x-integral: 

Lemma 27. There is a constant C 7 > 0 expressible in O, T and B such that 

' < C 7 S. 


R 7 ■■= At} - f V^{x)u'^{x) dx 

n = 1 «/ 


Proof. Since the product V^px is a smooth function on the domain we can invoke the mean- 
value theorem and find intermediate values 5^, such that 


rM 


6 S u(x^)p,(x^)- / U(Xl(C)p.(Xl(C)de 

fcelt ° 


The claim then follows by a change of variables. 


< 


5 ^ dx{VxPx)ixk){x^^i a^K-i) 


ken 


^ 5(6 - a) sup |I 4 (x)p 2 ,(x)|. 

£ceO 


□ 

































26 


HORST OSBERGER AND DANIEL MATTHES 


It remains to identify the integral expressions inside Ri to with those in the weak formu¬ 
lation (54). 

Lemma 28 . One has that 

rM 


22(0^€?A(ap'"oXl(OdC= i [ d4ulix)Yp"'ix)dx, 

^ Jfi 




R8-=tY, 


n=l 


rM 


^A{0iSizlf{0p"o^A{0d^- f {d.uir{x)p"{x)dx 

Jn 


(79) 

^ (80) 


for a constant Cg > 0 expressible in fl, T, B and 'Ey. 


Proof. The starting point is the relation (27) between the locally affine interpolants u\ and 
that is 

z2(6 = 21oXl(a (81) 

for all f e [0, M], Both sides of this equation are differentiable at almost every f e [0, M], with 

t5c22(0 = d,21oXl(ad5Xl(0- 


Substitute this expression for into the left-hand side of (79), and perform a change of 

variables x = X^(^) to obtain the integral on the right. 

Next observe that the x-integral in (80) can be written as 

r 1 

{d,u\f{x)p"{x)dx= / p" (82) 

JQ. Jo 

using (81). It hence remains to estimate the difference between the ^-integral in ([80|) and (82), 


respectively. To this end, observe that for each ^ e {fk,fk+^) some k e IJ, one has 
= V^fc+i and za(C) 6 [zk-ijZ^^i_]. Hence, for those 


1 - 


1 


1 - 




k— k 


If instead f e (^i,_i,^fe), then this estimate is satisfied with the roles of z7, i and z7_i inter- 

2 ^“r 2 ^2 

changed. Consequently, 

nM pM 1 

id^^Anozmp''oxi{f)df- ( 5 ^ 22 ) 2 (^)__^//oxi(c)d^ 


^BS Y, 4 

fcelt 


’ ~f£. ~f£. 

^fc+i ^k-i 


Zk+l 


- 1 


Zk-l 


de 

Zk-k 


- 1 


^/c+i 


which is again at least of order 0{6‘i), as we have seen before in (66). 


□ 


Proof of (62). Combining the discrete weak formulation (64), the change of variables formulae 
(79) and (80), and the definitions of i?i to Rs, it follows that 

Mx nM o rM 

e2.A ^BRs + BtY\ ^AiOS^zliOp'" o xm df+- zl{0i3^?^)^{0p" o XHf) df 

n=l'J 0 ^JO 

+ f v^(x) {uaIx (x)p\x) dx - 2 H” 

Jq 

7 7 

^bYR? ^bY 
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This implies the desired inequality (62). 


□ 


We are now going to finish the proof of this section’s main result, Proposition 18 


Proof of Proposition 18 Owing to (61) and (62), we know that 

rT 


10 


ri'{t) f p{x){uA}r{t,x)dx + r]{t)l- [ p'"{x)do,{{uA}l){t,x) + 3p"{x){doo {uA}^)‘^{t,x) dx 
Jq ^ Jn 


< 


Vj;{x) {uA}r {t,x)p'{x) dt 

2 

6i,A + 62,A ^ C(t + . 

To obtain (54) in the limit A ^ 0, we still need to show the convergence of the integrals to their 
respective limits, but this is no challenging task anymore: Note that ( |5^ implies 

dx {ua}x —> dxU:t: strongly in L^([0,r] x O), (83) 

hence {dx {ua}^.)^ converges to (dxU^)'^ in L^([0,T] x O). Furthermore, we have that 


dx{{uA}l) = ‘2{uA}rdx{uA}r -= dx(ul) 


(84) 


in L^([0,T] X O). Here we used (83) and that {ua}^. converges to uniformly on [0,T] x H 


due to (47). Hence, (83) and (84) suffice to pass to the limit in the second integral. Finally 


remember the weak convergence result in (46), {ua}^ —> m* in 7^2(f2) with respect to time, hence 
the convergence of the first and third integral is assured as well. 

□ 


5. Numerical results 

5.1. Non-uniform meshes. An equidistant mass grid — as used in the analysis above — leads 
to a good spatial resolution of regions where the value of is large, but provides a very poor 
resolution in regions where is small. Since we are interested in regions of low density, and 
especially in the evolution of supports, it is natural to use a non-equidistant mass grid with an 
adapted spatial resolution, like the one defined as follows: The mass discretization of [0, M] is 
determined by a vector S= (Co, 6, • ■ •,■Cjc-i, Ck), with 0 = Co < Ci < ’'' < Cif-i < = M 

and we introduce accordingly the distances (note the convention C-i = Cif-i = 0) 

(5k = Ck+1 - Ck-L, and 4 = ^ (<^ 10 +i +'5fc-i) 

for K e and k B respectively. The piecewise constant density function u e 
corresponding to a vector x e is now given by 

u(x) = Zk for a;„ 1 < x < x„ , 1 , with Zk = - - -• 

' ' (‘“2 «+ 2 ’ r 1 — T 1 

2 *^2 

The Wasserstein-like metric (and its corresponding norm) needs to be adapted as well: the scalar 
product (•, -fg is replaced by 

<v,w>^= 2 4'yfcWfc and llvll^ = <v, v>^. 
fc€l+ 

Hence the metric gradient Vy/(x) e of a function / : ^ K at x e yy is given by 
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Figure 1. Evolution of a discrete solution ua, evaluated at different times 
t = 0, 0.002, 0.012, 0.04 (from top left to bottom right) 


Otherwise, we proceed as before: the entropy is discretized by restriction, and the discretized 
information functional is the self-dissipation of the discretized entropy. Explicitly, the resulting 
fully discrete gradient flow equation attains the form 




— y""! 




(85) 


5.2. Implementation. Starting from the initial condition x^, the fully discrete solution is cal¬ 
culated inductively by solving the implicit Euler scheme (85) for x^, given In each time 


step, a damped Newton iteration is performed, with the solution from the previous time step as 
initial guess. 

5.3. Numerical experiments. In the following numerical experiments, we fix fl = (0,1). 

5.3.1. Evolution of discrete solutions. In a paper of Gruen and Beck [5], the authors analyzed, 
among other things, the behaviour of equation ([^ on the bounded domain (0,1) with Neumann¬ 
boundary conditions and the initial datum 


u^{x) = (x — 0.5)'^ -I- e, X e (0,1), with mass M = 0.0135, 


( 86 ) 


with e = 10“^. This case is interesting insofar as the observed film seems to rip at time t = 0.012. 
Figure shows the evolution of ua for K = 400 and t = 10“^ at times t = 0,0.0022, 0.012, 0.04, 
the associated particle flow is printed in figure [^left. 

5.3.2. Rate of convergence. For the analysis of the scheme’s convergence with initial datum 
with e = 10“^ from (86), we fix r = 10“^ and calculate solutions ua to our scheme with 


K = 25,50,100,200,400. A reference solution is obtained by solving (85) on a much finer 


grid, which is A = {K^J;Tref) with K^ef = 1600 and r = 5 • 10 In figurej^right, we plot the 
L^(r2), and L“(n)-norms of the differences |uA(t, •) ~ time t = 10“^. It is 


clearly seen that the errors decay with an almost perfect rate of 6^ccK 


5.3.3. Comparison with a standard numerical scheme. For an alternative verification of our 
scheme’s quality, we use a reference solution that ist calculated by means of a structurally 
different discretization of 0- Specifically, we employ a finite-difference approximation with step 
sizes Tref and h^^f = {b — a)/K,.e( in the t- and x-directions, respectively. More precisely, with 
tn ■= uTref, and with Xfc for fc = 0,..., AT^ef being the AT^ef + 1 equidistant grid points in [a, b], 
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Figure 2. Left: Associated particle flow of ua for initial datum (|86|). Right: 
Rate of convergence, using K = 25,50,100,200,400 and r = 10”^. 
are evaluated at time t = 10“^. 


The errors 




Figure 3. Left: Rate of convergence, using K = 25,50,100,200,400 and t = 
10“^. The discrete solutions are compared with a reference solution of the 
scheme in (87), and the errors are evaluated at time t = lO""^. Right: Loss of 


mass preservation using a standard finite-difference scheme and of (86) with 
e = 10“^, 10“^ and 10” 


- 1-5 


the numerical approximation ? 
— for given vector = (ug“^. 


u{tn\Xk) of 0 is obtained — inductively with respect to n 


) by solving the fully-implicit difference equation 


u — u 


n—1 


”^ref 


= U 


D 


ref ‘ 


+ dI 


ref ' 




(87) 


where u" = (mq ,..., and D)g£ are standard finite difference approximations of the ith 

derivative with equidistant steps /ij-ef- The product of two vectors in (87) shall be understood 
to act component-by-component. The boundary conditions 0 are enforced using values at 
“ghost points” in the obvious way, that is 


I L IL 1C It 1C 1C 1C 1C 

U-I—Uq, U_2 — Ui, + l UKrel+2 — 
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To produce reference solutions for the examples discussed below, the scheme above is imple¬ 
mented with itTref = 6400 spatial grid points and a time step r^ef = 5 • 10“®. At a given time 
T = NTref, the respective reference profile x Uref(T, cc) is defined via piecewise linear interpo¬ 
lation of the respective values 

Before comparing the advantages and disadvantages of our scheme and the reference scheme 

Using from (86) with £ = 10“^ instead 


5.3.2 


in (87), let us repeat the experiment of section 
of £ = 10“^, we plot the T^(U), A^(U), and L®(U)-norms of the differences \uA{t, •) — Uref(t, •)! 
at time t = 10“^ in figure [sj/lefb Obviously, the new experiment confirms the rate of convergence 
6^ccK~'^ gained in section [5.3.2| 

The following remarks concerning both schemes can now be made: 

• Computational cost: Using Newton’s method to solve both approximations (85) and (87), 

the finite-difference scheme is more efficient, unsurprisingly. The reason for this is the 
complicate structure of the Jacobian matrix of (x), whereas the Jacobian matrix of 

the right hand side of (87) is easy and quick to implement. Numerical experiments show 
that the finite-difference scheme can be approximately 5-times faster than our scheme, 
using the same values for K and r. 

• Conservation of mass: It is generally known that standard numerical schemes as the 

finite-difference approximation in (87) do not preserve mass. Depending on the initial 
datum, the loss or gain of mass can decrease quickely with time and causes inaccurate 
solutions. In figure 3 fright, we plot the relative change of mass Ure{(t, x) dx — I| 

with M = u^(x) da; for £ = 10“^, 10“^, 10“^ and t e [0,10“"*]. One can observe that 

the preservation of mass of solutions to the finite-difference scheme is seriously harmed 
in case of smaller choices of e. This is why we used e = 10“^ in the second experiment 
for the rate of convergence, since smaller values for e produce reference solutions whose 
change of mass yield to significant distortions of the L^-errors. 

• Conservation of positivity: In general, one can expect positivity of the discrete solution 
to the Hnite-difference scheme starting with a sufficiently positive initial function. This 
situation changes dramatically if one considers initial densities with regions of small 
values or even zero values. Take for example the initial datum in (86) with £ = 0. 
Then the solution to the scheme in (87) — again using = 3200 spatial grid points 
and a time step r^ef = 5 • 10“® — contains negative values after the very first time 
iteration and finally loses any physical meaning after some more iterations. In contrast, 
our scheme can still handle the case when e = 0 in (86), although one usually has to 
assume strict positivity for initial values in our approach. 

Conclusively, our scheme has a major advantage in comparison with standard numerical solvers if 
one is interested in a stable and structure-preserving discretization for ([^ . Moreover, the slightly 
plus of the finite-difference scheme and of similiar approximations discussed in the first point — 
less computational cost — is invalidated by the fact that one needs much finer discretization 
parameters compared to our structure-preserving scheme to gain solutions with an adequate 
physical meaning. 


Appendix A. Appendix 

Lemma 29 (Gargliardo-Nirenberg inequality). For each f e one has that 

ll/llcva(n) < (9/2)'/3||/||^?(a)ll/lli4"a)- (88) 

Proof. Assume first that / 5= 0. Then, for arbitrary x < y, the fundamental theorem of calculus 
and Holder’s inequality imply that 

< X / 1 •/(z)^''^l/'(2)|dz < - y\^^‘^\\f\\]!^a)\\f\\LHQ)- 

^ J X ^ 

Since / > 0, we can further estimate 

l/(^) - f{y)\ < ^ (3/2)2/3|:, - 
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This shows (881 for non-negative functions /. A general / can be written in the form / = /+ — /_, 
where /+ 5= 0. By the triangle inequality, and since ||/+||jLi(n) l/lHqn), 

ll/llci/6(n) ^ ll/+llci/6(n) + ll/-llci/6(n) < 2(3/2)2/3||/|| 

L2(n) ll/ll m{Q)- 

This proves the claim. 


Lemma 30. For each p ^ 1 and x G with z = Z 5 [x], one has that 

(x^^i - ^ {xk - xoY- 


2 (t)'- 


K6l 


1/2 


□ 


(89) 


Kel 


1/2 


Proof. The first equality is simply the definition 0 ofz„. Since trivially —ai^_i ^Xk — Xq 

1 /2 

for each k e , and since p — 1 > 0, it follows that 


EC 


X^,i — X, 


,_i)^ < {XK-Xof ^ 2 


.,_i) = {xK-Xo)P. 


□ 




1/2 


Kel 


1/2 
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